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ABSTRACT: The construction of a Supersymmetric Grand Unified Model based on the 
Pati-Salam gauge group is briefly reviewed and the low energy consequences of the derived 
asymptotic Yukawa quasi-unification conditions are examined. In the framework of the 
resulting Constrained Minimal Supersymmetric Standard Model, the cosmological relic 
density of the bino-like LSP is calculated and the results are explicitly compared with 
micrOMEGAs. In addition to the Cold Dark Matter constraint, restrictions on the parameter 
space, arising from the Higgs boson masses, the SUSY corrections to 6-quark mass, the 
muon anomalous magnetic moment and the inclusive decay b — > 57 are, also, investigated. 
For fi > 0, a wide and natural range of parameters is allowed. On the contrary, the ri < 
case not only is disfavored from the present experimental data on the muon anomalous 
magnetic moment, but also, it can be excluded from the combination of the Cold Dark 
Matter and BR(6 — ► 57) requirements. 
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1. Introduction 

This review is based on Refs. [1]. We briefly describe a Supersymmetric (SUSY) 
Grand Unified Theory (GUT) which predicts a set of asymptotic Yukawa quasi-unification 
conditions (YQUCs) in sec. 2 and then we show how we can restrict the parameter space 
of the resulting Constrained Minimal Supersymmetric Standard Model (CMSSM) using a 
number of Cosmo-Phenomenological requirements (sec. 3-4). Particular emphasis will be 
given in the Neutralino Relic Density calculation in sec. 3. Some conclusions and open 
issues will close this presentation in sec. 5. 

2. Yukawa Quasi-Unification 

The motivation of our construction is given in subsec. 2.1 and a brief description of 
the model building follows in subsec. 2.2. Finally, in subsec. 2.3, the mass parameters of 
the resulting CMSSM versions are derived. 

2.1 The Yukawa Unification Hypothesis 

We focus on a SUSY GUT based on the Pati Salam (PS) gauge group (SU(4) C ®SU(2) L 
®SU(2) R = 4 C 2 L 2 R ) described in detail in Ref. [2]. The relevant field content of the Model 
is presented in the Table 1, where the components, the representations, the transformations 
and the extra global charges of the various superfields are, also, shown. The left handed 
quarks and leptons superfields of every generation r are accommodated in a single pair 
of superfields F r ,Ff. The PS symmetry can be spontaneously broken down to Standard 
Model (SM) one through the vacuum expectation values (vevs) which the superfields H c , H c 
acquire in the direction of the right handed neutrinos, fjj^H- 
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In the simplest realization of this model, as it is proposed by Antoniadis and Leontaris 
in Ref. [3], the electroweak doublets H\,H 2 are exclusively contained in the bidoublet 
super field h, which can be written: 

h = ( h 2 h! ) , with hi = | ^- j and h 2 = [ j • (2-1) 

With these assumptions, the model predicts Yukawa unification (YU) at GUT scale, Mqut 
(Mqut is determined by the requirement of gauge coupling unification): 

h t (M GVT ) = h b (M GVT ) = h T {M GVT ) = y 33 , (2.2) 

since the Yukawa coupling terms of the resulting version of MSSM originate from a unique 
term of the underlying GUT, as follows: 

V33 F 3 h F 3 C 3 y 33 (Hjir 2 Q t c + hJit 2 Q b c + hJit 2 L r c ), since H 1[2] = h 1[2] . (2.3) 

Here Q,L 3 —ir 2 F 3 T are the well known SU{2)l doublets quarks and leptons of the third 
family and the self explanatory substitutions (t, b, r) in the position of (u 3 , d 3 , e 3 ) appeared 
in the Table 1 are made. Also, t 2 is the second Pauli matrix, 9 means "belongs to" and 
the brackets are used by applying disjunctive correspondence. 

However, applying this scheme in the context of the CMSSM [4] and given the top 
and tau experimental masses (which naturally restrict the tan/3 ~ 50), the pre-diction of 
the (asymptotic) YU leads to a contra-diction with the experimental bounds on 6-quark 
mass. This is due to the fact that in the large tan (5 regime, the tree level 6-quark mass, 
mj, receives sizeable (about 20%) SUSY corrections [5] Amb(MsusY), evaluated at a SUSY 
breaking scale Msusy (see sec. 2.3). These corrections arise from sbottom-gluino (mainly) 
and stop-chargino loops [5, 6, 7] and have the sign of /j, (with the standard sign convention 
of Ref. [8]). Their size is a ldoop exact quantity, as it is explained in Ref. [7]. Hence, 
for \i > [<] 0, the SUSY corrections drive the corrected (which is essentially the so called 
running [7]) 6-quark mass, m£, at a low scale Mz, 

m c b {M z ) = m b (M z ) (l + Am 6 (M SU SY)) (2.4) 

a lot above [a little below] than its 95% confidence level (c.l.) experimental range. This is 
derived by appropriately [1] evolving the corresponding range for the pole 6-quark mass, 
mb(mt,) up to Mz scale with a s (Mz) — 0.1185, in accord with the analysis in Ref. [9]: 

m b (m b ) = 4.25 ± 0.3 GeV m b (M z ) = 2.88 ± 0.2 GeV (2.5) 

Concluding, for both signs of fi the YU assumption leads to an unacceptable 6-quark mass. 
The discrepancy is less sizable for /j, < 0, so we can say that YU "prefers" [10] /i < 0. 

The usual strategy to accommodate the latter discrepancy is the introduction of several 
kinds of nonuniversalities in the scalar [11, 12, 13] and/or gaugino [10] sector of MSSM, 
without an exact conservation of the YU. On the contrary, in Ref. [1], this problem is 
addressed in the context of the PS GUT, without need of invoking the departure from the 
CMSSM [4] universality. The Higgs sector of the model is extended so, that the electroweak 
Higgs in Eq. (2.3) receive subdominant contributions from other representations, too. As 
a consequence, a moderate violation of the YU is obtained, which can allow an acceptable 
6-quark mass, even with universal boundary conditions. 
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Table 1: Relevant Field Content 
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2.2 Model Construction 

From the discussion of the previous section, we can induce that a small deviation from 
the YU will be enough for an appropriate prediction of 6-quark mass when \x < 0, while 
for n > 0, a more pronounced deviation is needed. According to the first paper in Ref. [1], 
this can be achieved by introducing some extra higgs fields, as follows: 

Since the matter fields product has the following structure under 4 c 2l2^: 

F 3 F 3 C B (4, 2, 1) © (4, 1, 2) = (1, 2, 2) © (15, 2, 2) (2.6) 

and 155(/(4) B lsu(3)> it is possible the addition of two higgs superfields: 

h', h! e (15,2,2). (2.7) 

The field h' can couple to F3F3 and h! can give mass to the color non singlets through 
a mixing term mh'h' with m ~ Mgut — 2 x 10 16 GeV, in accord with the imposed 
global symmetries (see Table 1). Other possible mixing term suppressed by the string scale 
M s ~ 5 x 10 17 GeV, being non renormalizable, is: 

H c H c h'h/M s 3 (l5 ©15, 1©1, (1© 3) ®(10 3)), (2.8) 

with the latter structure under 4,^21,2^, since for the 2 participants of the product, we get: 

H C H C [h'h] G (l5 © 1 [15], 1 [1 © 3], 1 © 3) . 

In the term of Eq. (2.8), there are 2 couplings as regards the SU(2)r: A singlet, X\{H C H C )\ 
h'h (since 1(g) 1 = 1), and a triplet, X^(H c H c )^h'h (since 3©3 = 10305). As it turns out, 
the singlet coupling provides us with an adequate deviation from the YU for /j, < 0. The 
necessary deviation for \i > can be obtained by a further enlargement of the Higgs sector, 
such that contributions arising from renormalizable terms are allowed. By introducing two 
additional Higgs fields: 

(15,1,3) (2.9) 

an unsuppressed coupling X' 3 <fih'h (again since 3©3 = 1©3©5) can be constructed. This 
overshadows the coupling X^(H c H c )^h'h. <p is introduced to give superheavy masses to the 
color non singlets in <j) through a term <f><j). Summarizing, the important mixing terms are: 

mh'h 1 , \ 1 [ 3 ](H c H c ) 1 { 3 }h'h/Ms and \' 3 4>h'h (Xi^, X' 3 : dimensionless constants) . 

Using the transformation properties of the several fields indicated in the Table 1, and 
taking into account the unitarity of the SU (2) l r groups, we conclude that the previous 
mixing terms correspond to the following invariant quantities under 4 c 2^2 j r, respectively: 

mTr (h'iT 2 h' T iT 2 ^ , (h'iT 2 (H cT H cT ) 1[3] h T iT 2 ^ and A' 3 Tr (h' ir 2 (j)h T iT 2 ^j , (2.10) 

where the traces are taken with respect to the SU(4) C and SU(2)l indices. 

At GUT scale, during the spontaneous breaking of PS gauge symmetry, the fields 
H c , H c and <p acquire vevs in the SM singlet direction. Therefore, 

[(H cT )(H cT )) m = v 2 Hc (t H c, 1, ^) and (<P) = V<p (V 5 , 1, , (2.11) 
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where VH c ,<f> ~ AfcuT and the three-fold structure of these vevs as regards the A c 2l2r 
group is ordered in the parenthesis by commas, with 

T H c = diag(o, 0,0,l), T 15 = ^=diag(l,l,l,-3) and a 0[3] = diag(l, [-]l) . (2.12) 

Expanding the superfields in Eq. (2.7) as linear combination of the 15 generators T a of 
SU(4:) C with the normalization Tr(T a T b ) = 5 ab and denoting the SM singlet components 
with the superfield symbol, the following identities can be easily proved: 

Tr (h!iT 2 h! T iT 2 ^ = h?ir 2 ti 2 + h?ir 2 h 2 + 0{T a T a , a < 15), (2.13) 
Tr (h'i T2 (rg c] , 1, <7 3[0] ) h T iT 2 ) = (h?ir 2 h 2 - [+]h[iT 2 h! 2 ) [-V3/2] , (2.14) 

where the notation of Eqs. (2.11) was applied. Inserting Eqs. (2.11) in Eqs. (2.10) and 
employing Eqs. (2.13)-(2.14), we get the mass terms: 

mh?iT2 (ti 2 + a 2 h 2 ) + m (h? + a^J^j ir 2 h' 2 , (2.15) 
where the the mixing effects are included in the following coefficients: 
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(2.16) 



It is obvious from Eq. (2.15) that we obtain 2 combinations of superheavy massive fields: 

h 1 A- n h ■ 

h[, H[ and h' 2 , H 2 , where H[ = 1 = , i = l,2. (2.17) 



The electroweak doublets, remaining massless at GUT scale, can be interpreted as "or- 
thogonal" to H[ direction: 

Hj = Z$MM (2.18) 



Solving Eqs. (2.18) and (2.17) as respect hi and h[, we obtain: 

-atHi + H' 
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^1 + |ai| 2 ^1 + H 2 

Consequently, the Yukawa couplings terms, in contrast to Eq. (2.3), take the form: 

y 33 ^3 h Fl + 2y' 33 F 3 h' F§ 3 y 3 3 ^3(^2 + 2pT l5 h> 2 h, + 2pT 15 /i' 1 ) FS o 3 
1 pOL 2 j ■s/S j . 1 - pai/y/S T . c 1 + v^pai T . c 



(2.19) 



2/33 



\ \A + M 2 " z " y/l + |o!i| 2 " _1 " ' x/TTlail 
where p = 2/33/2/33 with < p < 1 and only color singlets components of h\ are shown. The 
fields H!- being superheavy mass states, contribute neither to the running of renormalization 
group equations (RGEs) nor to the masses of fermions. Consequently, the asymptotic exact 
YU in Eq. (2.2) can be replaced by a set of asymptotic YQUCs: 



ht(M GTJT ) : h b (M GUT ) : h T (M GUT ) 



1 - pa 2 /V~3 



1 — pa\j\f?> 



1 + \[?>pa\ 



\A + l"i| 
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V! + l«2| 2 

The deviation from the YU can be estimated, by defining the following relative splittings: 
Sh b[T] = {h b[T] (M GVT ) - h t (M GVT )) /ht(M GVT ). (2.21) 
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FIGURE 1: A Yukawa couplings running from M GVT to M z for m c b (M z ) = 2.68 GeV, /j, > [<] 0, 
tolsp = 297 [372] GeV and A? 2 = 1.5 [0] corresponding to tan/3 ~ 59 [50] (left [right] plot). 



2.3 CMSSM With Yukawa Quasi-Unification 

Below the Mqut, our particle content reduces to this of MSSM. We assume universal 
soft SUSY breaking scalar masses, too, gaugino masses, M 1 / 2 and trilinear scalar couplings, 
A at M G ut- Therefore, the resulting MSSM is the so called CMSSM [4] supplemented 
by a suitable YQUC from the set in Eq. (2.20). With these initial conditions, we run the 
MSSM RGEs [14] between M GUT and a common SUSY threshold M SUSY ^ (to^to^) 1 / 2 
(ti 5 2 are the stop mass eigenstates) determined in consistency with the SUSY spectrum. At 
Msusy we impose radiative electroweak symmetry breaking, evaluate the SUSY spectrum 
and incorporate the SUSY corrections to b and r masses [6, 7, 11]. The latter (almost 4%) 
lead [14] to a small de[in]-crease of tan/3 for [i > [<] 0. From Msusy to Mz, the running 
of gauge and Yukawa couplings is continued using the SM RGEs. 

For presentation purposes, M 1 / 2 and too can be replaced [14] by the lightest SUSY 
particle (LSP) mass, tolsp and the relative mass splitting between it and the lightest stau 
T2, Af 2 . For simplicity, we restrict this presentation to the Aq = case (for Aq ^ see 
Refs. [1, 15]). Our parameter "transmutation" is shown schematically, as follows: 

sign/x, M 1/2 , to , A) sign/x, to LS p, A f2 (A f . 2 = (m f2 -to LS p)/to LS p) 

We restrict ourselves to monoparametric YQUCs derived from the general Eq. (2.20) 
with real 0:1,2 ■ It is convenient to define the quantity: 

c = pai/s/3 with — 1/3 < c < 1 

For any given to^(M^) in Eq. (2.5) with fixed top, mt(mt) = 166 GeV and tau, m c T {Mz) = 
1.746 GeV masses, we can determine the parameters c and tan/3 at Msusy so, that the 
YQUC in Eq. (2.20) is satisfied. However, for the plots in Figs. 1, 2 and 3 we fix 
ml(Mz) = 2.68 GeV, since this value turns out to overlap the allowed areas for fi > or 
give the less restrictive upper bound for \i < on the tolsp — Af 2 plane (Figs. 4 and 5). 
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FIGURE 2: The mass parameters m a, mo, M 1 / 2 and Msusy versus m^sp form%(Mz) — 2.68 GeV, 
\i > [<] and A f2 =1.5 [0] corresponding to tan/3 ~ 59 [50] (left [right] plot). 

The YQUC predicting appropriate values for m b (Mz) relates the minimal Higgs con- 
tent of the model described in sec. 2.2 to the sign of fi, as follows: 

• For fi > the fields of Eq. (2.9) are needed. Viable CMSSM spectrum is obtained 
for ol\ = —(X2 ~ 1 in Eq. (2.16). Therefore, Eq. (2.20) implies: 

ht(M GVT ) : h b (M GVT ) : h T (M GVT ) = |1 + c| : |1 - c\ : |1 + 3c|. (2.22) 

In the left panel of Fig. 1, we present a typical Yukawa couplings running from Mgut to 
Mz for mLSP = 297 GeV and A T2 = 1.5. In the left panel of Fig. 2, the values of m^, 
mo, Mi/2 and Msusy versus mLSP are, also, displayed for the same A T2 . We observe that 
M 1/2 < m and 2m LSP ~ m A for m LSP < 350 GeV. For 2.68 GeV < m c b (M z ) < 3.09 GeV, 
the parameters c, 5hu T ] and tan (3 range, correspondingly: 

0-19 > c > 0.13, 0.32 > Shr = -5h b = 2c/ (1 + c) > 0.23 and 60 > tan/3 > 57. 

• For fi < 0, the inclusion of the fields in Eq. (2.9) can be avoided. Interesting CMSSM 
spectrum is obtained for a\ = «2 ~ 0.1 in Eq. (2.16). Therefore, Eq. (2.20) implies: 

h t (M Gl]T ) : h b (M GVT ) : /i r (M GUT ) = |1 - c| : |1 - c\ : |1 + 3c|. (2.23) 

In the right panel of Fig. 1, we present a typical Yukawa couplings running from Mqut 
to Mz for ttilsp = 372 GeV and Af 2 = 0. In the right panel of Fig. 2, the values of m-A, 
mo, Mi/2 an d Msusy versus mLSP are, also, displayed for the same Af 2 . We observe that 
M 1 / 2 > m and 2m LS p > m A . For 2.68 GeV < m c b (M z ) < 3.09 GeV, the parameters c, 
<5/i6[ T i and tan (5 range, correspondingly: 

0-01 < -c < 0.06, 0.04 < -Shr = -4c/(l - c) < 0.25, 5h b = and 50 > tan/3 > 47. 

Worth mentioning is finally, that in our models (in contrast to the "traditional" 
CMSSM version [4]) tan/3 is not a free parameter but a prediction of the applied YQUC. 
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3. Neutralino Relic Density 

In the context of CMSSM, the LSP can be the lightest neutralino. It naturally arises 
as a Cold Dark Matter (CDM) candidate [16]. We require its relic density, f^LSP^ 2 ) n °t to 
exceed the upper bound derived from DASI on the CDM abundance at 95% c.l. [17]: 

Ocdm/i 2 < 0.22 (3.1) 

For both signs of fx an upper bound on tolsp can be derived from this requirement. Some 
generic features of JIlsp^ 2 calculation in the CMSSM are exploited in subsec. 3.1, and 
particular applications to the model under consideration are exhibited in subsec. 3.2, 3.3. 

3.1 General Considerations 

In most of the CMSSM parameter space, the LSP is almost a pure bino and Q^sph 2 
increases with ttt-lsp- Therefore, Eq. (3.1) sets an upper approximate limit on its mass: 
m-LSP ^ 200 GeV. However, as it is pointed out in Refs. [18, 19], a substantial reduction of 
r^LSP^- 2 can be achieved in some regions of the parameter space, thanks to two reduction 
"procedures": The Apole effect (APE) and the coannihilation mechanism (CAM). The 
first is activated for tan (3 > 40 [~ (30 — 35)] with [i > [<] 0, where the presence of a 
resonance (2mLSP — fnA) in the A Higgs mediated LSP' s annihilation channel is possible. 
On the other hand, CAM is applicable for any tan/3, both signs of fi but it requires a mass 
proximity between LSP and the next-to-LSP, which turns out to be the lightest stau, T2, 
for tan/3 > 10 [20, 21] and not too large values for A [22] or m [23, 24]. 

Our model gives us the opportunity to discuss the operation of both reduction "pro- 
cedures". As it is induced from Fig. 2, for fi > there is a significant region with APE, 
while for fi < the CAM is the only available reduction "procedure". Furthermore, for 
fi < 0, the CAM is more strengthened, since due to the larger mLSP, more coannihilation 
channels are kinematically allowed than for \x > 0. For this reason some technical details 
for each reduction "procedure" will be presented separately in the two following subsections 
for fj, > and /i < respectively. 

We calculate r^sp^- 2 ) using micrOMEGAs [25], which is one of the most complete publicly 
available codes. This includes accurately thermally averaged exact tree- level cross sections 
of all possible (co)annihilation processes, and loop QCD corrections to the Higgs couplings 
into fermions. The results of this code are checked in Refs. [1, 15] with another public 
package DarkSUSY [26] (not the newest version [24]) appropriately combined [1] with the 
code used in Ref. [14]. We found good agreement when the Higgs couplings are treated 
at tree level. This agreement persists, even with loop QCD corrections to these couplings 
[29], provided an artificial m^rrib) is used in the defaults of DarkSUSY, in order to mimic 
these corrections. 

In this talk, a new comparison is presented with an improved version of the code (let 
name it, GLP) presented in Ref. [14]. A first, model independent improvement concerns 
the freeze out procedure [27]. This is renewed, using variable values for the number of 
relativistic degrees of freedom, g*, whose a precise estimation is obtained, employing the 
tables included in micrOMEGAs package. Details on other specific improvements together 
with the comparisons will be displayed in the next subsections. 
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3.2 Including the A-Pole Effect 



Table 2: Comparison of Olsp^ 2 Calculations 





Model Parameters 






(ji > 0, m c b (M z ) = 2.68 GeV, A = GeV) 




tan P 


59 


58.9 


58. 


8 1 


M 1/2 (GeV) 


712 


663 


616 


m (GeV) 


1298 


1484 


1638 


m A (GeV) 


560 


510 


462 


m LS p (GeV) 


319 


297 


276 


Af 2 


1.0 


1.5 


2.0 


^Lsph 2 Calculations 


Code : 


micr— 
□MEGAs 


GLP 


micr- 
OMEGAs 


GLP 


micr- 
OMEGAs 


GLP 


r^ cc (GeV) 


66 


64.4 


59 


58.4 


54 


52.6 


r3 CD (GeV) 


34 


32.5 


31 


29.8 


29 


27.0 


tsusy (GeV) 




24.9 




23.3 




21.5 


f^LSP^ 2 tree 


0.124 


0.117 


0.123 


0.116 


0.122 


0.115 


Qhsph 2 QCD 


0.221 


0.218 


0.223 


0.221 


0.222 


0.221 


^LSP^ SUSY 




0.270 




0.271 




0.267 



In the absence of the CAM, ^lsp^ 2 is inverse proportional to the thermally averaged 
product of relative velocity times the x ~ X cross section, (va^x)- ^ n the presence of the 
APE, this is enhanced, due to the s-channel A exchange, s(A), to the down type fermion- 
antifermion pairs. 

For a reliable calculation of Olsp^ 2 in this regime, two points have to be taken into 
account: First, since the low velocity expansion of (va^x) breaks down [27] in the vicinity 
of poles, the full phase-space integration for the s(A[_ff]) channels with fermions to final 
states has to be performed, for masses in the interval [20] 0.65 < 2m^/mA\H] < 2 (the H 
contribution is p-wave suppressed [20, 28]). Second, a careful treatment of the relevant 
A- fermions couplings gAff j and the corresponding decay widths is, also, indispensable [25]. 
Since these are proportional to the corresponding fermion mass, m/, a rather accurate 
estimation of their tree level [29], g^f f {M), loop QCD corrected [29], g^f \M) and SUSY 
resummed [7] QCD corrected g A l j i j Y (M) value at a scale M can be achieved, if rrif is 
replaced by an effective fermion mass. More explicitly, 

{g'XffiM) , m f = m f (m f ) 
g^fiM) , m f = m/(M)(l + Am f (M SVSY )) A^ 2 D (M) (3.2) 
ff |USY (M)) m/ = m/(M )Aj/ 2 D (M) 



917lf t&n(3= { 



2M W 

with Aqcd = 1 [given in Ref. [29] (applying the itia 3> m/ limit 1 )] for / = r [b] (g is the 



1 We are grateful to Dr. P. Ullio for bringing to our attention this point. 
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FIGURE 3: n LSP h 2 versus m hS p for m c b (M z ) = 2.68 GeV, fj, > [<] and A f2 = 1.5 [0] (left [right] 
plot). The thick (faint) solid line is obtained by micrOMEGAs with (without) loop QCD corrections 
[co annihilations] while the thick (faint) crosses by our code. The stars are obtained by our code, 
including SUSY corrections, too. The upper bound on J7lsp^ 2 (= 0.22) is also depicted. 

SU(2)l gauge coupling). Inserting these effective couplings calculated at ijia in the tree 
level formula for the ^4-decay width [29], we reach the micrOMEGAs results in the case of 
tree level, T^ ee and QCD corrected, r^ CD A-dec&y width, within a 5% accuracy as it is 
shown in the Table 2. An estimation, also, can be done for the SUSY improved value of 
the A-decay width, r^ USY (not included in the current version of micrOMEGAs). 

For the inputs of Table 2, the relative contributions of the annihilation processes to 
the Q^sph 2 calculation are as follows: 

XX ~~ ^ bb [rf] 93 [6]% at tree level, 

with a 6% de[in]-crease at loop QCD level and a 8% de[in]-crease at SUSY improved loop 
QCD level. In the same Table, we compare, also, the corresponding values for the tree 
level, ^LSP^ 2 |tree and loop QCD corrected, ^lsp^ 2 |qcd> ^lsp^ 2 with micrOMEGAs. An 
agreement within a (4-0)% accuracy is achieved. Note that the various couplings in our 
code are calculated at 2mLSP scale. Moreover, the 45% increase of Q^sph 2 because of 
the loop QCD corrections to <M//'s and ^4-decay width (firstly noticed in Ref. [25]) is 
impressively reproduced. An estimation can be done, as well, for the result on Olsp^ 2 , 
if SUSY corrections are included, ^lsp^ 2 |susy- An almost 18% increase is expected. 
This statement can not be reliably checked through micrOMEGAs, since the user is able to 
introduce these corrections only to ^4-decay width and not to gAff 's [25], too. 

A further comparison is displayed in Fig. 3 (left plot), where we depict fiLSP^ 2 versus 
w-LSP for fj, > 0, Af 2 = 1.5 and mb(Mz) = 2.68 GeV. fiLSP^ 2 is calculated by micrOMEGAs 
[our code] with (thick solid line [crosses]) or without (faint solid line [crosses]) the inclusion 
of the loop QCD corrections to gAff couplings. The stars are obtained by our code, 
including the SUSY corrections, too. 
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3.3 Including Bino-Stau Coannihilations 



Table 3: Included Feynman Diagrams 



Initial State 


Final State 


Interaction Channels 


XX 


// 

J J 


s(h),s(H),s(A),s(Z) 




( f: Fermions) 


t( f 1 9), u( f 1 9) 






s(h),s(H),s(A),t(xt),u(xf) 




hA,ZH 


s(A),a(Z),*(v?),«(x?) 


YT9 
A ' ^ 


rh, tH, tZ 


s(t), t(fi 2 ),u(y?) 




tA 


V / ; v v J- / 5 \A^l J 




TJ 


s{t)Mt 2 ) 




u T H~ ,u T W~ 


s{r),t(xf),u(i) T ) 


t 2 t 2 


TT 


f (*?),«(*?) 


r 2 r 2 


hh, hH, HH, ZZ 


s(h),s(H),t(f lt2 ), U (f lt2 ),Pl 




AA 


s(fe),a(fr),*(fi),«(fi),PI 




H+H-,W + W~ 


a(/i),s(fO,s(7),a(Z),*(i/ T ),PI 




H±WT 


s(/0,s(#),*(*t) 




AZ 


s(/0,a(#),*(fi),u(n) 




11, lZ 


i(f 2 ),«(f 2 ),PI 




ti,bb 


s (/i), a (fr), a (7),s(z) 




TT 


s(h),s(H),s( 7 ),s(Z),t(x) 



Bino-stau coannihilations come into play, when Af 2 < 0.25. Q.j J sph 2 is not any more in- 
verse proportional to (va^) but to an effective (va) which includes, in addition, (wfaxlfa]) 
with a weight factor exp (— pJA-^Xp 1 ), where Xp 1 = mLSp/^F ~ 25 and Tf the freeze out 
temperature [27]. Consequently, for given m^sp, A^ 2 regulates the degeneracy amount. 
The strongest possible reduction is achieved for Af 2 = 0. 

The computation of f^LSP^ 2 m this regime, can be realized by using exclusively the 
low velocity expansion of (va). The thermal average has been performed, following the 
Ref. [21]. The included set of (co) annihilation processes (in accord with the tables in Refs. 
[28]) is shown in the Table 3 (we apply the notation of Ref. [14] and PI stands for "point 
interaction"). The relevant matrix elements are evaluated with the help of the FeynCalc 
package [30]. Our results are compared with micrOMEGAs, in the Table 4 for 3 test points. 
Besides the values of ^lsp^ 2 and Xp 1 , we display, also, the contribution of all the channels, 
beyond 0.3%. As we can observe, our differences are small in most cases. 

The importance of CAM, in deriving an upper bound on tolsp from Eq. (3.1), can 
be easily concluded from Fig. 3 (right plot), where we depict ^lsp^ 2 versus ttt-lsp for 
|U < 0, Af 2 = and m^(Mz) = 2.68 GeV. Olsp^ 2 is calculated with (thick solid line 
[crosses]) or without (faint solid line [crosses]) the inclusion of bino-stau coannihilations 
using micrOMEGAs [our code]. We see that the reduction of S7lsp^ 2 caused by the CAM is 
dramatic (a factor of 10). This increases the upper bound on tolsp derived from Eq. (3.1) 
from about 200 GeV to 900 GeV. 
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Table 4: Comparison of ft^sph 2 Calculations 





Model Parameters 








(jjt < 0, rrf b {M z ) = 2.68 GeV, A = GeV) 






tan (3 


51.2 


50.0 


49.4 




M 1/2 (GeV) 


1877 


859 


596 




m (GeV) 


957 


463 


357 




m A (GeV) 


466 


214 


130 




m LS p (GeV) 


865 


386 


264 






0.00 


0.05 


0.10 




^lsp^ 2 Calculations 


Code : 


micr— 
OMEGAs 


GLP 


micr- 
OMEGAs 


GLP 


micr- 
□MEGAs 


GLP 


Processes Which Contribute More Than 0.3% 


Process 


Contribution (%) 


XX // 


0.4 


0.5 


8.4 


10.6 


25.7 




30.5 


xx - 


0.8 


0.7 


9.4 


10.1 


18.6 




19.6 


XX^HZ 


0.4 


0.2 


5.1 


3.8 


11.2 




10.4 


xx ^ 


0.4 


0.3 


5.2 


4.3 


11.0 




9.3 












0.3 




0.3 


XX ^ 










0.3 




0.2 


Xf 2 -> r/i 


1.4 


0.4 


8.0 


6.9 


7.1 




6.1 


XT2 — > ri7 


3.0 


2.8 


4.8 


4.6 


2.6 




2.4 


Xf 2 -> 


3.0 


3.3 


4.7 


5.3 


2.7 




2.8 


XT 2 -> f T #~ 


5.8 


6.0 


8.5 


9.1 


4.4 




4.4 


XT 2 -»■ U T W~ 


0.6 


0.5 


1.5 


1.3 


0.9 




0.7 




7.9 


7.7 


11.9 


12.3 


6.4 




6.3 


Xf 2 -> rZ 


2.7 


2.7 


4.8 


5.0 


2.8 




2.7 


^2T2 — > TT 


19.0 


19.2 


5.5 


5.7 


0.9 


0.8 


T2f| — > /i/l 


8.0 


8.2 


4.2 


3.9 


0.8 




0.6 


f 2 f 2 * -»• HH 


3.0 


3.2 


0.9 


0.9 








f 2 f 2 * -► 


3.0 


3.2 


0.9 


0.9 








f 2 f 2 * -► ff+fl- 


5.6 


6.2 


1.3 


1.5 








f 2 f 2 * -► ff^W* 


0.4 


0.4 












f 2 f* _^ vf+vf- 


13.8 


14.4 


4.4 


4.6 


0.6 




0.6 


f 2 f| -► zz 


7.3 


7.7 


2.3 


2.4 


0.3 




0.3 


T2T2 -> 77 


5.7 


5.3 


1.7 


1.7 


0.3 




0.2 


f 2 f| -► 7Z 


3.0 


3.2 


0.7 


0.8 








f 2 f 2 — > 


4.1 


3.0 


4.5 


3.5 


0.9 




0.6 




25.97 


26.03 


25.51 


25.48 


24.82 




24.84 


^LSP/1 2 


0.221 


0.205 


0.221 


0.221 


0.220 




0.216 
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4. Phenomenological Constraints 

A lower bound on iblsp can be derived by imposing a number of phenomenological 
constraints [10, 13, 20, 23]. These result from: 

i. The Higgs boson masses. The relevant for our analysis is the 95% c.l. LEP bound 
[31] on the lightest CP-even neutral Higgs boson, h mass 

m h > 114.4 GeV (4.1) 

which gives a [almost always the absolute] lower bound on ttilsp for /j, < [>] 0. The SUSY 
contributions to are calculated at two-loop by using the FeynHiggsFast [32] program 
included in micrOMEGAs package [25]. 

ii. The deviation 5a^ of the muon anomalous magnetic moment measured value, 
a^, from its predicted in the SM, a^ M . The latter is not yet stabilized mainly due to 
the instability of the hadronic vacuum polarization contribution. According to the most 
updated [35] evaluation of this contribution, the findings based on e + e~ data and on r- 
decay data are inconsistent with each other. Combining these results with the recent 
experimental measurements on [33], we obtain the following 95% c.l. ranges: 

a) 11.3 x 10" 10 < 5a^ and b) 5a^ < 56.1 x 10~ 10 e+e^-based (4.2) 
a) - 11.6 x 10~ 10 < 5a^ and b) 5a^ < 30.4 x 10~ 10 r-based (4.3) 

The SUSY contribution to the 5a^ is calculated by using the formulae of Ref. [34] in accord 
with micrOMEGAs [25] and the result is posi[nega]-tive for /i > [<] 0. A lower bound on 
mLSP can be derived for //>[<] from Eq. (4.2b) [(4.3a)] and an optimistic upper bound 
for fi > from Eq. (4.2a) which, however is not imposed as an absolute constraint. 

iii. The inclusive branching ratio of b — ► s 7 , BR(6 — ► s 7 ). Taking into account the 
recent experimental results [36] on this ratio and combining appropriately the experimental 
and theoretical involved errors [1], we obtain the following 95% c.l. range: 

a) 1.9 x 10~ 4 < BR(6 -► s 7 ) and b) BR(6 -► s 7 ) < 4.6 x 10" 4 (4.4) 

To calculate BR(6 — > 57) we used an updated version of the code contained in the current 
version of micrOMEGAs [37]. This code represents a complete update respect the one used 
in the first paper of Ref. [1]. The SM contribution is calculated using the Refs. [38]. The 
charged Higgs boson [SUSY] contribution is evaluated by including the next-to-leading 
order QCD [SUSY resummed] corrections and tan/3 enhanced contributions from Refs. 
[39]. A lower bound on m^sp can be derived for fj, > [<] from Eq. (4.4a) [(4.4b)] with 
the latter being much more restrictive. 

5. Conclusions-Open Issues 

Allowing ml(Mz) to vary in its 95% c.l. range, Eq. (2.5) and using the previous 
Cosmo-Phenomenological restrictions, we can delineate the parameter space of the model 
on the mLSP — Af 2 plane in Figs. 4 and 5. For simplicity, we do not show bounds from 
less restrictive constraints, Eq. (4.2b) and (4.4a) [(4.1)] for \i > [<] 0. 
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™ LSP (GeV) m LSp (GeV) 



FIGURE 4: The allowed area on the TOLSP-Af 2 FIGURE 5: Disconnected areas on the m^sp- 

plane for fj, > and 2.68 GeV < m c b (M z ) < A f2 pJane for fi < and 2.68 GeV < mg(M z ) < 

3.09 GeV. Shown is, also, the [optimistic] upper 3.09 GeV allowed by various restrictions indi- 

bound from Eq. [(4.2a)] (3.1) with 2.61 GeV < cated in the plot. 
m b (M z ) < 3.15 GeV, [dashed] dotted line. 

For fi > 0, the allowed space of parameters (depicted in Fig. 4) is wide as regards both 
Af 2 (because of the APE) and m-LSP (because of the CAM) range. Especially, 

(From Eq. (4.1) : ) 159 GeV < m LSP < 720 [371] GeV ( : From Eq. (3.1) [(4.2a)]) 

(taking the less restrictive values of the various constraints). The changes on the upper 
limit of the allowed area from a possible inclusion of the SUSY corrections (sec. 3.2) in 
the Q,LSph 2 calculation are almost invisible. On the contrary, let a s (Mz) vary in its 95% 
c.l. experimental range 0.1145 — 0.1225 with the corresponding range of m c h (Mz) to be 
(2.62 — 3.16) GeV, the allowed area can further enlarge until the dotted line in the Fig. 4. 

On the other hand, for \x < 0, the competition of the various Cosmo-Phenomenological 
constraints is more dangerous and, finally, disastrous for this case, since: 

(From Eq. (4.3a) : ) 387 GeV < m LS p < 873 GeV ( : From Eq. (3.1)) 

and m LS p > 1305 GeV ( : From Eq. (4.4b)) 

Despite the strong presence of the CAM (Af 2 < 0.05) we are, evidently, left without 
simultaneously allowed region, as it is shown in Fig. 5. 

Constructions similar to this presented in sec. 2.2 may be useful for other SUSY GUTs, 
too. The ,50(10) or SU(5) SUSY GUTs in their simple realization lead to complete [12, 13] 
or b — t [10, 40] YU, with the latter being viable only for fi < and consequently, quite 
disfavored from the bounds of Eq. (4.2). Also, the large values of tan/3 predicted in these 
models, can enhance the neutralino detection rates with universal [41] or non universal 
[10] asymptotic gaugino masses. Furthermore, the extra higgs fields used in sec. 2.2 have 
interesting consequences to the inflation mechanism, as it is pointed out in Ref. [42]. 
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